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We  prove  consistency,  stability  and  convergence  of  the  point  vortex  approximation  to  the.  2-D 
incompressible  Euler  equations  with  smooth  solutions. "Wa  first  show  thatfthe  discretization  error  is 
second-order  accurate.  Than  show  thar’the  method  is  stable  in  ^  norm.  Consequently  the  method 
converges  in  ^  norm  for  all  time.  The  convergence  is  also  illustrated  by  a  numerical  experiment.  j<j2jp  y ,  (  . 


We  prove  consistency,  stability,  and  convergence  of  the  point  vortex  approxi¬ 
mation  to  the  incompressible  Euler  equations  in  two  dimensions.  The  vortex 
method  for  this  problem  (see  e.g.  Chorin  [5])  has  two  parameters,  the  grid  length, 
h,  and  the  smoothing  length  (or  “blob  size”),  5.  In  computational  practice  people 
usually  take  8  -  h  or  5  =  0  (see  [13])  but  the  convergence  theorem  of  Hald  [10] 
and  subsequent  work  needed  the  hypothesis  that  8  »  h,  leaving  these  cases  of 
computational  practice  in  doubt.  In  this  paper  we  prove  convergence  with  8  =  0 
under  hypotheses  that  are  otherwise  similar  to  those  of  Hald.  This  convergence  is 
also  illustrated  by  a  numerical  experiment.  Of  course,  we  can  not  claim  that  the 
point  vortex  method  is  better  than  the  vortex  blob  method  or  vice  versa;  our 
purpose  is  to  show  that  the  point  vortex  method  works  in  principle.  In  future 
work  [9],  we  shall  treat  the  case  8*0  but  8  not  much  larger  than  h.  We  shall 
also  treat  some  of  the  proposed  three-dimensional  vortex  methods;  see  [11]. 

The  proof  follows  that  of  Hald  with  an  improvement  (lp  instead  of  l2  norms) 
introduced  by  Beale  and  Majda  [3].  We  show  that  the  method  is  formally 
second-order  accurate  by  showing  that  certain  sums  are  second-order  accurate 
approximations  to  the  corresponding  Biot  Savart  integrals.  Actually,  we  get  an 
asymptotic  expansion  for  the  truncation  error  that  would  allow  us  to  use  Strang’s 
device  in  [15]  and  prove  convergence  for  low-order  accurate  methods  in  high 
dimensions.  That  is  deferred  to  a  later  paper.  We  then  prove  lp  stability  for 
1  <  p  <  co  by  using  Lp  bounds  for  a  related  continuous  singular  integral 
operator,  much  as  was  done  by  Hald  and  by  Beale  and  Majda.  Except  for  this  Lp 
bound,  our  proof  is  elementary  and  selfcontained. 
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In  vorticity  formulation,  the  two-dimensional  in  viscid  incompressible  Euler 
equations  are  (see  [10]) 

(1)  o>,  +  (u  *  v)w  =  0, 

(2)  S’  «(*>')  =  (K{x -y)u(y,t)  dy , 

(3)  *(*)  ■ 

It  is  known  (see  e.g.  [16])  that  for  the  initial  vorticity  m(jc,0)  =  u°(x)  which  is 
smooth  and  localized  («°  e  S,  the  Schwartz  class;  see  [8])  a  solution  of  (1),  (2) 
exists  for  all  times  («(•,  t)  e  S  for  all  t  >  0).  The  Lagrangian  (or  “flow  map”) 
coordinate,  £,  is  defined  implicitly  by  the  particle  trajectories,  y{£,  /),  which 
satisfy 

(4)  ^yU,t)  =  u(y(i,t),t),  y(£,  0)=£. 

If  w°  e  S,  then  this  change  of  variables  is  smooth.  Since  the  flow  is  incompress¬ 
ible,  det (dy/d£)  =  1  and  therefore  the  mapping  £  -*  y  has  a  smooth  inverse 
£  =  £(y,  /).  The  point  vortex  method  approximation  is  to  replace  (2)  and  (1)  by 

(5)  Vj(t)  **  A2  £  K(xj(t)  -  xk(t))«*» 

k+j 

and 

(6)  Xj(t)  *  Vj(t), 

where  xy(0)  -  h  -  j,  j  =  (jv  j2)  eZJ,  h  €  R  is  the  grid  length,  and  «y  *= 
<o°(xy(0)).  A  consequence  of  localized  vorticity  is  that  particles  far  away  move 
very  little.  In  particular,  for  any  T  we  have 

lim  sup  |j>(€,/)  -  I!  “  0. 

«|ir, 

>ST 


From  (1)  it  follows  that 


«(M€.  0,  t)  «  «(£,0)  -  «°({). 


The  exact  particle  trajectories  yXt)  are  the  Lagrangian  images  of  the  original 
grid;  y/t)  -  y(£Jt  t).  We  write  £}  - h-j  so  that  w(yy(0. 0  -  «°(£y)  -  «y-  The 
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discrete  lp  norm  is  defined  by 

ll/ll '  -  A2£l/,l'- 

j 

Recall  that  v(t)  is  the  approximate  velocity  field  defined  by  (5).  We  also  write 
Uj  =  u(yj(t),  t ),  the  exact  solution.  Our  result  is 

Theorem  A.  For  all  T  and  4  <  p  <  oo,  there  is  a  C{T,  p)  such  that 

(7)  | \x(t)-y(t)HC(T,p)h\ 

(8)  \\v(t)  -u(t)\\zC(T,p)h2. 

Theorem  B.  Suppose  u°  e  S,  4  <  p  <  oo,  omf  T  >  0  are  given,  and  that 

(9)  ||*(/)-t(0IU*2‘1/4  f<>r  Os*  £7. 

Then  there  exists  an  A(T,  p)  such  that 

(10)  HO  -  «*(/)  II  £  A(T,  p)||*(0  -  y(t)  II  +  A(T,  p)h\ 

for  0  g  t  £  T,  where  A(T,  p)  is  a  nondecreasing  function  in  T. 

Remarks.  1.  The  assumption  «°  e  S  can  be  replaced  by  a  finite  degree  of 
differentiability  and  decay. 

2.  The  h2~1/A  in  (9)  can  be  replaced  by  ||x(/)  -  y(t)  ||/<o  g  <>(/r)A  for  any 
positive  function  <j>(A)  with  ^(A)  -»  0  as  A  -*  0. 

3.  By  a  longer  consistency  argument  that  uses  the  asymptotic  error  expansion 
(17)  in  the  manner  of  Strang  [IS],  the  restriction  4  <  p  <  oo  can  be  replaced  by 
1  <  p  <  oo. 

4.  Our  Lemma  1  to  follow  is  in  the  spirit  of  the  remark  of  Beale  [2]  on  the 
consistency  of  a  3-D  vortex  blob  method. 

5.  The  K*9iniqiu»s  in  this  paper  also  apply  to  the  deterministic  vortex  method 
of  Cottet  and  Mas-Gallic  [6]  for  the  incompressible  Navier-Stokes  equations. 

6.  Our  result  does  not  apply  to  the  ill-posed  vortex  sheet  problem  where 
analytidty  seems  to  be  necessary  for  convergence;  see  [4],  . 

We  choose  to  present  a  proof  that  does  not  depend  on  smoothing  the  kernel 
AT  as  in  the  blob  method.  We  mention,  however,  that  due  to  the  incompressibility 
and  smoothness  of  the  flow  ™*p  y({),  the  point  vortex  method  can  be  identified 
with  the  blob  method  with  blobs  of  size  0(A);  see  [11].  Consequently,  the  blob 
t^hniqiwn  can  be  more  or  less  applied  directly  resulting  in  a  shorter  but  less 
self  contained  proof.  We  do  not  present  this  method  to  emphasize  that  no 
additional  smoothing  is  necessary. 
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The  statement  of  direct  practical  interest  is  Theorem  A,  but  it  is  more 
convenient  to  prove  Theorem  B.  Let  us  indicate  how  Theorem  A  follows  from 
Theorem  B. 

Proof  of  Theorem  A  from  Theorem  B:  For  a  given  h  and  p,  define  T*  by 

T*  =  inf{t:  ||x(f)  -y(OII  *h2~t/4}. 

It  suffices  to  show  that  if  h  <;  ^  exp{  —4 A(T,  p)  •  T),  then  T*  >  T.  For  in  that 
case  ( T*  >  T),  we  have  |j jc( / )  —  y(t)  ||^  h2~1/4  for  all  t  £  T.  It  follows  from 

(10)  that,  for  t  <;  T, 

a  l*(«)  -  >(«)  I  S IM<>  -  «(')  II  S  A(T,  P)lx(,)  -y(l)l  +A(T,  p)h\ 

which,  via  Gronwall’s  inequality,  proves  (7)  with  C  =  eA'T.  We  get  (8)  with 
C  =  AeA'T  from  (7)  and  (10).  Now,  if  7"*  s  T,  we  have,  for  all  /  g  T*, 

^  MO  -HOI  ^IM 0  -  “(Oil  £  A(T*,p)||x(r)  -  y(0  II  +  A(T,p)h2, 

by  (10),  which  implies  that  MO  -  y(t)  ||g  exp {A(T*,  p)T*}h 2  ior  t  <,T*  by 
Gronwall’s  inequality.  But  A(T*,  p)  is  a  nondecreasing  function  of  T*.  By  the 
choice  of  h,  we  conclude  that 

MO  -HOI  £  exp{A(T,p)T}h2g  \h2~ V* 

for  all  t  g  T*,  which  contradicts  the  definition  of  T*. 

Proof  of  Theorem  B:  First  we  show  that  the  y}  are  close  to  satisfying  (5)  and 
(6).  More  precisely,  define  the  residual,  or  “truncation  error”,  p,  by 

(11)  Pj(t)  -  u(yj(t),  t)  -  h2  £  K(yj(t)  -  yk(t))uk. 

k+j 

We  claim: 

Lemma  1  (Consistency).  There  is  a  constant,  D(T),  such  that  ||p(0  II £ 
D(T)h2  for  all  t  £  T. 

Proof  of  Lemma  1:  In  (2)  we  change  from  y  to  £  as  variable  of  integration 
and  set  x  -  y/t).  Using  incompressibility,  we  have 


«U<0.0-  jK(yAt)-y(e,t))u°(e)d(. 
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Thus  p  is  the  error  in  trapezoid  rule  integration.  Therefore  the  lemma  would  be 
routine  but  for  the  singularity  of  K.  To  handle  this  we  split  p  into  its  near  field 
part,  a,  and  its  far  field  part,  t.  For  this,  we  use  a  smooth  cut-off  function,  /(s), 
of  s  e  R,  satisfying  f(s)  *  1  if  s  g  1  and  f(s)  *  0  if  s  £  2.  Let  g(s )  =  1  —  f(s) 
be  the  complement  of  /.  We  split  up  pj  into  Oj  +  rJt  where  0  <  q  <  1  and 


(12) 


and 


« >-  /*U-y(0)«#U)/(*'%-*M 

-  E  *(yj  -  yk)<*°(tk)f(*-%  -  «*0 

k*j 


(13) 


tj  -  fK(yj -y(S))o°(t)g(h-%  -  £|)  di 
-hlY,K(yj-  yk)»0Uk)g(h-%  -  «*l). 

k+j 


For  notational  convenience,  we  temporarily  take  j  =>  (0,0)  so  that  £  *  (0,0) 
and  we  write  y'  for  y'(0),  the  derivative  of  y  with  respect  to  £,  etc. 

To  estimate  the  near  field,  we  must  characterize  the  near  field  behavior  of 
tf(y(0)  -  Xl))«°(0  for  small  £ 


Lemma  2.  For  any  N,  we  have 

(14)  *(->-(£))  «°(£)  -  ^  +  «0(f )  +  *i(«)  +  •••  +0(\t\N), 

where  b  is  a  2x2  matrix  and  m„  is  homogeneous  of  degree  n  in  (:  m„(a£)  « 
a"in„({)  for  all  a  e  R. 

The  denominator  does  not  vanish  since  det (y')  -  1.  The  terms  b  and  m„ 
depend  on  t,  but  they  as  well  as  the  error  bound  implied  in  (14)  are  uniform  in 
t  £  T.  Since  the  m„  are  actually  rational  functions  of  £  a  need  not  be  positive. 
In  particular,  m„(£)  is  an  even  or  odd  function  of  £  depending  on  whether  n  is 
even  or  odd. 


Proof  of  Lemma  2:  Note  that  since  y(£  t)  and  w°(£)  are  C°°  functions  of  £ 
we  have  the  expansions 

(15)  y(£) -m+si  +  r(u)  +  —  +  0(1*1*) 

and 

(16)  «°(£)  -  «°(0)  +  «*'£  +  «°"(£  £)  +  •••  +  0(|£|") 


\ 
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for  any  N.  The  numerator  of  K  has  a  simple  expansion,  so  we  need  only  study 
the  denominator 

1 _ 1 

|y(f)2l  rt  +  yl 

But,  using  (15),  we  get 


2(y1'  ‘  i)  •  V(j,  t)  +  2  (y2' - l)  () 

\y'i\2 


where  the  successive  terms  have  the  orders  and  smoothness  claimed.  Once  we 
have  an  expansion  for  K  using  (15),  we  can  multiply  it  by  (16)  to  get  the  desired 
expansion  for  Ku,  which  proves  the  lemma. 

Consider  the  integral  term  in  (12).  Since  /(|{|)  is  even,  we  have 


and 


/»»«(€)/(*  *lflW€  =  0  for  nodd. 


Also,  the  substitution  ij  =  h~9£  gives,  for  n  even, 

/"*„(£)/(  A-* lil)  d(  -  a„A<"+2>*  where  an  *  di\. 

A  similar  analysis  holds  for  the  discrete  summation  terms  of  (12).  The  odd 
order 'terms  vanish.  For  the  even  terms  we  have  a  little  trick.  The  sums  may  be 
written  as 

*#0. 

where 

SM  =  E«,(i)/(A‘lt|). 

k*  0 

Now,  since  S„  is  a  sum  with  finitely  many  non-zero  terms, 

4rS,(h)  *(1-9)  I  «I,(*)V/(*1-*|*0  *  *'«[*! 


-  (1  -  £  m(,(A1-»k)/(A‘-^|)(A1-*)2, 

*1*0 
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where  f(y)  =  V^/flyD  •  \y\.  The  sum  is  a  trapezoid  rule  approximation  to 

K  =  fm*(y)f(y)  dy- 

The  integrand  is  a  C®  function  of  y  since  V/  vanishes  at  the  origin  where  m,  is 
singular  and  vf  has  compact  support.  It  is  well  known  (e.g.  [1],  [7],  p.  300)  that 
the  trapezoid  rule  for  such  integrals  has  an  infinite  order  of  accuracy,  that  is 

(A1-*)2!  =bn  +  0(hN)  forany  N. 


This  shows  that 


±Sn(h)  =  (1  -  q)h-^-^bH  +  0(h»). 

Using  this  and  the  integral  representation 

S„{h)  -  SB(1)  -  dh 

we  obtain  .. 

S„(/r)  =cn  +  bHh-<”+2*l-‘»+  0(hN)  forany  N, 

where  c„  is  a  constant  and  bK  —  -b^n  +  2). 

Putting  all  this  together  and  returning  to  the  general  case  j  *  0  gives 

°j  -  0*2  +  cA(yj,  t)h 4  +  •  •  • } 

+  (d2(y,,  r)A**  +  dA(yj,t) h*9  +•••)  +  0(hN)  forany  N. 

The  coefficients  c„  and  dn  are  functions  of  derivatives  of  u  and  y  at  the  point  y} 
which  decay  rapidly  at  infinity  if  a  e  S.  Comparing  the  expansion  corresponding 
to  |  <  q  <  1  with  the  expansion  corresponding  to  0  <  q  <  we  conclude  that 
d„  -  0  for  n  even.  That  is, 

Oj-  c2{yj,t)h2  +  c4(yJtt)h4  +  •••  +0(hN)  forany  N. 

•  / 

The  theory  of  numerical  integration  tells  us  that  t}  -  0(hN)  for  any  N ,  since 
the  integrand  is  without  singularities  (thanks  to  g)  and  decays  at  infinity  (thanks 
to  <•>).  Combining  this  with  our  estimate  for  or,,  we  conclude  that 

,  •' 

(17)  Pji'h) -ca((y,f)AJ  +  c4Uy,/)*4  +  •••+©( hN)  forany  N. 
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Since  the  coefficients  cH(y)  are  rapidly  decaying  as  y  -»  oo,  summing  over  j 
shows  that  ||p||  is  0(h2)\  this  proves  Lemma  1. 

We  note  that  it  is  possible  to  obtain  second-order  accuracy  by  using  the 
oddness  of  the  most  singular  part  of  the  kernel  and  using  straightforward 
estimates  involving  the  midpoint  rule  (see  [11]). 

We  turn  now  to  proving  (10).  We  use  the  following  estimate:  For  any  mesh 
function  fp 

(18)  A2l/,l'£A2El/*l'“  ll/ll'; 

k 

thus 

i/,1  3  A-2/'ll/l|. 

In  particular,  applying  (18)  to  /  -  x(t)  -  y(t)  and  using  (9)  gives  us 

(19)  |*,(<)  - •>>(<)  | 

We  conclude  that,  for  p  >  4,  the  left  side  of  (19),  the  error  in  any  particular 
particle  position,  is  small  relative  to  A,  the  typical  distance  between  particles. 

Another  useful  fact  is  Young’s  inequality  (e.g.  [8],  p.  13).  If  A  *  {aJk}  is  an 
N  x  N  matrix,  and 

Mr  «  max£|<j.A|  and  Mc  =•  max£|ay*|, 

J  lc  k  j 

(the  subscripts  r  and  c  stand  for  row  and  column,  respectively)  then 

Mil  £  max(Afr,  Me). 

Similarly,  if  L(x,  y)  is  an  integral  kernel  with 

»up  f\L(x,y)\dy  and  Me  -  sup  f\L(x,  y) \dx, 

x  "y  y  *  x 


then 

£  max(Af,,  Me). 


We  are  now  ready  to  estimate  v  -  u.  Using  (S)  and  (11),  we  have 
py(0  -“/(*)  £(*(•*>“**)  ~  K{yj-yit))uk  +  Pj- 

k+j 
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Using  a  Taylor  expansion  for  each  component  of  K,  we  can  write  this  as 
Vj{t)  -Uj(t)  =  A2£  VK{yj-yk)((xj-yj)  ~(xk-yk))uk 

k+j 

+  A2  £  V2K(zj  -  zk)\Xj  -yj~(xk-  yk)]osk  +  Pj, 

k+j 

where  V2K(x)[y]  stands  for  the  matrix  of  second  partials  acting  as  a  quadratic 
form  on  the  argument  that  appears  in  square  braces.  From  the  mean  value 
theorem,  Zj  is  somewhere  on  the  line  segment  connecting  Xj  to  yjy  and  zk  is  on 
the  segment  connecting  xk  to  yk.  (This  is  an  abuse  of  notation  since  Zj  depends 
on  k  and  vice  versa.)  It  follows  from  (19)  that 

(20)  Zj=yj  +  o(h). 

Decompose  v  -  u  as  a  sum  of  Tx  +  T2  +  T3  +  p,  where 

Tx  =  h2l  £  vK(yj-  yk)a,\(xj-  yj), 

' k+j  ' 

Tz  =  h2  £  vK(yj  -  yk)uk(xk  -  yk), 

k+j 


T3  =  h2  £  v 2K{zj  -  Zk) [xj  -yj-{xk-  yt)] 

k+j 

Lemma  1  gives  an  lp  bound  for  p.  We  shall  derive  lp  bounds  for  each  of  the  Tt  of 
the  form  ||7)||  g  C\\x  -  y||  for  /  =  1, 2, 3. 

Using  the  definition  (3)  of  K  and  (20),  we  bound  the  second  derivative  term 
by 


V2K(zj-zk)\  g 


I zj  z*'3  \yt  -yk  +  °(h)\3 


Since  £  is  a  smooth  function  of  y,  and  since  \y}  -  yfc|  £  c(/)|£;  -  £*|  ^  c(t)h  for 
j  *  k,  the  right  side  is  bounded  by  C(f)/(|£y  -  £*|3).  Estimating  the  sum  by  the 
corresponding  integral  gives 


k+j 
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where  a(t)  >  0.  With  this,  we  have,  using  (19)  and  Young’s  inequality. 

Fall  £  A(p,  T)hl-'W\\x(t)  -  y(t )  ||. 

To  handle  Tx  and  T2  we  use  Hald’s  trick  of  dividing  space  into  Lagrangian 
boxes.  These  are  defined  by  (here  y  =  (yt,  y2)  and  j  -  (jv  j2)): 

B}(t)  =  {ytt,  0  :  nvaxdir  -  hjx\,  ||2  -  hj2\)  g  \h) 


and 


Bj(t)  =  {y(4,f)  :  max(|{x  -  hj2\,\i2  -  hj2 1)  g  i*}. 


The  smaller  boxes,  Bjt  seem  necessary  to  keep  points  in  one  small  box  separated 
from  points  in  neighboring  boxes.  Because  of  the  volume  preserving  property  of 

Hy), 


(21) 


(  dy  as  h2  and  f  dy  =  \h2. 

JBjO )  JBj(t) 


We  first  use  the  decomposition  into  boxes  to  bound  the  sums  in  T2: 

dj  =  h2X,VK(yj~yk)uk, 
k+j 

which  are  approximations  to 

/vAT(y,-y(0)«0(€)^. 

Since  Tx  is  a  diagonal  matrix  with  diagonal  entries  dj  acting  on  Xj  -  yp  the  lp 
bound  on  Tx  follows  from  the  boundedness  of  the  individual  dj.  First,  by  the 
smoothness  of  u  and  y,  and  the  approximate  oddness  of  K  (as  in  Lemma  1),  we 
have 


■if 

J 

A 


f  9  fC(yj-y)u(y,t)  dy 


(22) 


«°(*/)/  9*K(yj-y)4 

JBj(t )  n 

K(yj-y)4yi 


+  o(h) 


+  0(h)  -  0(h). 


y- 

1l7. 


'..r'V 
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Now,  we  integrate  over  each  box  (ignoring,  as  we  may  by  (22),  Bf)  and  sum  over 
the  boxes  to  get 


A  E  VK{yj-  yk)u>\  -  fvK(yj  -  y(())a°(()  d( 

k 

k*j 


*hj 


”  &j\ 


by  the  midpoint  rule  (see  [7]),  where  we  have  used  the  fact  that  \y{ij)  —  y(£)l  ^ 
a(t)|«y  -  (|  £  «{t)h  for  y(t)  *  Bj(t). 

Thus,  we  have 


\dj\=\h2'LvK{yJ-yk)ok  £C{t). 

1  k+j 


This  settles  Tt  but  leaves  T2.  For  this,  we  need  an  lp  bound  for  the  matrix  whose 
elements  are  hWK(yj  -  yk)uk.  But  multiplication  by  uk  is  bounded,  so  we  only 
need: 


Lemma  3  (Stability).  The  matrix  M  whose  elements  are 


ih^vKiyj- yk)  if  j  +  h, 

10 


has  the  lf  bound  ||iWf ||  C(f,  p). 

Proof  of  Lemma  3:  Suppose  /  is  any  grid  vector  and  g  *  Mf.  We  want  to 
show  that  ||g||  ^  C(t,  p)\\f\\.  For  this  (following  Hald),  we  introduce  functions 
of  y.  First,  F(y )  i^  defined  in  terms  of  /  by  F(y)  -  Afj  if  y  e  BJt  and  by 
F(y)  -  0  if  y  £  U jBj.  The  factor  4  is  such  that,  in  view  of  (21),  ||/||,  -  ||F|| ,  . 
In  the  same  way,  G(y)  is  defined  in  terms  of  g.  Define  a  piecewise  constant 
kernel 

Lh(y,y')  m  l*K(yj-yk)  if  y  e  Bj(t),  y’  e  Bk(t)  and  j  #  k, 

\  0  otherwise, 

and  the  corresponding  “exact”  kernel 


L{y,/)  -  VK(y-y'). 
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The  characteristic  function  for  Bj  is 


xj(y)  = 


if  y  e  Bp 

otherwise, 


and  the  characteristic  function  for  UjBj  is  x(y)  =  Y,Xj(y)- 

j 


xOO- 


fi 

io 


if  y  e  u Bj, 

j 

otherwise. 


Schematically,  the  proof  is  as  follows: 


ll*lli,  =  IK?II l,  =  4||x  *  L„  •  xF||L/  £  4C\\LF\\Lf  g  4C'||F||^  =  4C'||/||V 


where  LF(y)  =  j L(y,  y')F{y')  dy'.  We  must  establish  the  two  questioned  in¬ 
equalities.  The  second  of  them  is  the  Calderon  Zygmund  inequality  (e.g.  [11], 
p.  29). 

The  remaining  inequality  is  proved  by  showing  that 
||x(L*-L)x||t,£C(/). 


We  divide  the  kernel 

R(y>  y  )  =  x(y)(Lh(y<  /)  -  Hy,  y'))x(y') 

into  two  parts:  R  =  +  R2.  First  we  take  the  case  when  y  and  y'  are  in  the 

same  small  box: 

Rx(  y<  y')  =  /  R(y’  y’)  if  y*=Bj  and  ?  e  \  for  some  j, 

1 0  otherwise 

If  /?x(y,  y')  *  0,  then  Lh(y ,  y')  =  0;  thus 

£  E  JjLXjF\'dy  <;  £||L(x,F)II£, 

j  bj  j 

S  m\l£\\XjF\\lr  z  \\L\\l'\\F\\lp. 
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By  the  Calderon-Zygmund  inequality,  we  get  ||#{|||£  ^  B(p)  <  oo.  Now  we  look 
at  the  part  coming  from  different  boxes:  ' 

Rz(y,y')  =  x(^)(vA:(y>  -yk)  -  VK(y  -y'))x(y') 

if  ye  Bj  and  y'  e  Bk  and  j  #  k. 

For  R2,  we  use  Young’s  inequality.  The  usual  mean  value  argument  gives 

I vK(yj  -  yk)  -  VK(y  -y')l  S  (l>)  ~y\  +  I yk  ~  /l)  *  I ~  zk)\, 

where,  as  always,  yy  -  y  =  0(h)  and  yk  —  y’  =  O(h).  Now,  since  the  boxes  Bj 
and  Bk  are  separated  by  0(h),  we  have  the  bound 


/ 1 R2(  y>  y')  I  dy'sChZf  IV  2K(  Zj  -  zk)  \  dy' 


Ch  f 
J\i- 


{'ia«(Q*li  ~  i'l3 


di'SC(t). 


This  argument  is  symmetric  in  y  and  y'  and  so  it  also  gives  the  other  bond 
needed  for  Young’s  inequality, 


/l*2U  y')\dy$C{t). 

This  completes  the  proof  of  the  lemma  and  of  the  theorem. 

A  numerical  experiment  also  demonstrates  the  second-order  convergence  of 
the  point  vortex  method.  We  started  with  initial  conditions  (see  [12]) 

uo(x)  /(1-W2)7,  0$  |*|  *1, 

\0,  |*|  >1, 

and  solved  the  equations  (6)  using  fourth  order  Runge  Kutta  (see  e.g.  [7],  p.  346) 
with  h  =  At.  It  is  easy  to  compute  the  exact  velocity  field  and  Lagrangian 
coordinate.  The  velocity  field  does  not  change  with  time.  Particles  near  the  origin 
complete  one  rotation  at  time  t  =  4w,  while  the  particles  on  |x|  =  1  complete  one 
rotation  at  t  —  32ir.  The  velocity  error  is  E(h,  t)  =  ||i>  —  u||.  Figures  1  and  2 
show  the  velocity  error  in  maximum  and  /2  norms,  respectively,  for  h  = 
0.2,0.1,0.05.  The  order  of  accuracy  is  the  power,  r,  in  E(h,t )  ~  C(t)hr.  In 
Figure  3,  we  plot  r(h,  t)  =  \og2(E(2h,  t)/E(h,  t)\  which  is  an  estimate  of  r. 
For  sufficiently  short  times,  the  second-order  accuracy  (r  =  2)  is  borne  out.  In 
Figure  4  we  compare  E(h,  t)  to  ||p(/)||(see  (11)).  The  agreement  of  these  two 
indicates  that  the  method  is  stable. 


\ 


Figure  2.  /,  error*  in  velocity,  h  :  (-J0.2,( —  )0.1,  ( . . .  )0.05 


> 
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